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Abstract 

An experiment is proposed to excite the 'forbidden' ls2s^Si - ls2s^So magnetic dipole (Ml) 
transition at 1.557/xm in a collimated and slow atomic beam of metastable helium atoms. It is 
demonstrated that an excitation rate of 5000 s~^ can be realised with the beam of a 2 W narrowband 
d ' telecom fiber laser intersecting the atomic beam perpendicularly. A Doppler-limited sub-MHz 

CJ ' spectroscopic linewidth is anticipated. Doppler-free excitation of 2% of trapped and cooled atoms 
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J>^, may be realised in a one-dimensional optical lattice geometry, using the 2 W laser both for trapping 

and spectroscopy. The very small (8 Hz) natural linewidth of this transition presents an opportunity 

for accurate tests of atomic structure calculations of the helium atom. A measurement of the '^He 
> ■ 

. - He isotope shift allows for accurate determination of the difference in nuclear charge radius of 

Q both isotopes. 
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INTRODUCTION 



Measurements of level energies of low-lying states in helium provide very sensitive tests 
of basic theory of atomic structure. Although helium is a two-electron system, energy levels 
of the nonrelativistic helium atom can be calculated with a precision that is, for all practical 
purposes, as good as for nonrelativistic hydrogen. Relativistic corrections to these energies 
can be calculated in a power series of the fine structure constant a and have been calculated 
up to 0[o?\ Effects of the finite nuclear mass can be included in a power series of the mass 
ratio /i/M, where /i is the reduced electron mass and M the nuclear mass. QED corrections, 
both hydrogenic and electron-electron terms, up to O(a^) have been calculated as well. 
Effects of the finite nuclear size can be incorporated straightforwardly [1|, |2| . Present day 
theory aims at calculating higher-order corrections (in a, /i/M) and cross terms (such as 
relativistic recoil). The most difficult terms to date are the relativistic and QED terms of 
O(a^) and higher. Higher-order corrections are largest for low-lying S-states and therefore 
the most sensitive tests of atomic structure calculations can be performed for the 1 ^Sq 
ground state and the metastable states 2^So (lifetime 20 ms) and 2^Si (lifetime 8000s). 

Present-day laser spectroscopy on the low-lying S-states has several disadvantages. To ex- 
tract an experimental value for the ionisation energy of these states the transition frequency 
to a high-lying state has to be accurately measured and one has to rely on theoretical values 
of the ionisation energy of the upper state in the transition. Here the lifetime of the upper 
state and line shifts due to stray electric fields and laser power are limiting factors. For the 

t round state, excitation is difficult: one photon of 58 nm is required to excite the 2 ^Pi state 
I or two photons of 120 nm to excite the 2 ^Sn state Ml. For the 2 ^Sn and 2 ^Si metastable 
...,Q,.,aH-..j 

states. Two-photon spectroscopy is applied to excite the 3^D state j9[ and n^D states |10|. 

The most accurate transition frequency measurement to date is for the 2^Si 2^Po 
transition j^, with an absolute accuracy of 4 kHz. This measurement, however, does not 
provide an accurate measurement of the 2 ^Si ionisation energy as the 2 '^Pq ionisation energy 
is not known well enough. Moreover, recent measurements of the fine structure splitting of 
the 2'^P state have shown that experiment and theory do not agree at the 10 kHz level 11 1. 
This very recent finding is considered an outstanding problem of bound state QED and asks 
for independent measurements on other transitions with similar accuracy. Experimental 
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lan present-day theory for the 2 ^Sq and 2 ^Si ionisation energy, which is accurate 
and 1 MHz respectively. The theoretical accuracies represent the estimated 



accuracies for the ionisation energy of the 2 ^Sq and 2 ^Si states, deduced from measurements 
to highly excited states and relying on the theoretical calculations of the ionisation energies 
of these states, are 150 kHz and 60 kHz respectively. These values agree with but are more 
accurate t^ 
to 5 MHz 

magnitude of uncalculated higher-order terms in QED calculations. 

QED shifts largely cancel when identical transitions in different isotopes are studied. 
Measuring the isotope shift, the main theoretical inaccuracy is in the difference in the rms 
charge radius of the nuclei jl^ . Q| . As the charge radius of the ^He nucleus is the most 
accurately known of all nuclei (including the proton) isotope shifts measure the charge radius 



of the other isotope involved. In this way accurate determination of the 2 ^Si - 2 
isotope shift have allowed accurate measurements of the charge radius of ^He 3, ll^ and 
the unstable isotope ^He |l7|, challenging nuclear physics calculations. It has to be noted, 
lowever, that the measurement of the ^He nuclear radius has sofar not been reproduced 
3|- These measurements therefore primarily measure differences in charge radius. 

Elaborating on an idea of Baklanov and Denisov we propose direct laser excitation of 
the 2'^Si 2^So transition in a slow atomic beam or in an optical lattice. This transition 
has the advantage of an intrinsically narrow natural linewidth of 8 Hz and a wavelength of 
1.557 /im. Also, as the transition connects two states of the same ls2s configuration, the 
theoretical error in the transition frequency may be smaller than the error in the ionisation 
energy of the individual states The main disadvantage is that the transition is extremely 
weak; the Einstein A-coefficient for this magnetic dipole (Ml) transition is 14 orders of 
magnitude smaller than for the electric dipole (El) 2'^Si — 2'^P transition. In this paper 
we show that with 2 W of a narrowband fiber laser at 1.557yum we can excite more than 1 
in 10^ atoms in an atomic beam experiment or more than 1% of atoms trapped in a one- 
dimensional standing light wave. Present-day sources of metastable helium atoms [2^ easily 
provide sufficient atoms in either an atomic beam or in a trap to observe the transition. 



transition 



EXPERIMENTAL FEASIBILITY 

The magnetic dipole transition between the metastable 2'^Si state and the 1 ^Sq ground 
state determines the ~ 8000 s lifetime of the 2 ^Si state. The experimental and theoretical 
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FIG. 1: Relevant levels, transitions and decay rates (left) and simplified model (right) 



values of the rate constant for this transition, 1.10(33) x 10^^ s~^ [2l| and 1.272 x 10^'^s^^ 
13,01 respectively, agree very well. For the electronically similar transition from 2^Si to 
2^So, rate constants of 1.5 x 10^"^ s~^ [24] and 6.1 x 10^^ s~^ have been published. In 
this letter we will use the value 9.1 x 10~^ s~^, obtained by Pachucki applying the same 
formalism as used for the calculation of the 2'^Si to 1 ^Sq transition rate [23[. This value is 
assumed to be accurate at the 1% level. 

In order to evaluate the feasibility of detecting the transition, we need to choose an 

n 

experimental configuration. Several approaches can be considered |2^: spectroscopy in a 
discharge cell, on a thermal or laser-cooled beam, on a cold cloud in a magnetic or optical 
trap, and ultimately spectroscopy on atoms in an optical lattice. Here, we first consider a 
relatively simple but promising setup: spectroscopy on a laser-cooled and coUimated beam 
as produced on a daily basis in groups working on BEC in metastable helium |2o| . 

BLOCK EQUATIONS MODEL 

In order to describe the excitation of the 2'^Si 2^So magnetic dipole transition, we 
use a simplified set of optical Bloch equations (OBE's). The decay rate of the upper level 



is fully dominated by the two-photon (2E1) decay to the ground state: r2 = 51 s~^ \24\. 
Fig. n summarizes the relevant decay rates. The driven Rabi frequency of the transition 
(with atomic frequency LJ21) is denoted by Q, the detuning by A = ujiaser — '^21- 
We model this transition by the following set of OBE's: 



Pii = ^2iP22 - Tipii + i^(p2i - P12) (1) 



P22 = -r2P22 - ^^(P21 - P12) (2) 
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P'l2 = - (^^-^ + i^^ Pl2 + i^{p22 - pu) (3) 
P2i = -y ^ P2i-i^{p22- Pn) (4) 

As the ground state is not included, this set of equations does not have a steady state 
state solution except the trivial one (pn = p22 = Pu = P21 = 0). 

When we simplify this set of equations by neglecting A21 and Fi (which is certainly valid 
for times t <^ l/Fi ~ 8000 s), it can be solved analytically. Starting at t = with all atoms 
in the 2 ^Si state (pn = 1), and denoting r2 by T from now on, the result for the population 
P22 of the 2 ^Sq state is: 



P22(t) 



29? 



cosh I jl[r]2_f]2 + ^] t 



COSM/^[^^2 + f^2_^]^ 



(5) 



where fi^ = /(^ + ^)2 + A2 - ^2 + ^2 and 9^ = 9'' + A' 



In this Letter, two limiting cases are studied, both valid for short time (t ^ 2/r = 40 ms): 
the weak excitation limit {1/t ^ Q ^ r/2) and the strong excitation limit {Q ^ 1/t ^ 
r/2). For the weak excitation limit, Eq. reduces to: 

"sin(A t/2) 



P22(t)^^sin2(A t/2) 



A t/2 



(6) 



In the strong excitation limit, we can derive a simple expression for the upper level 
population time-averaged over the oscillations of the cosine in Eq. El 



(P22) 



2(^2 + A2) 



(7) 



BROADENING EFFECTS 



Eqs. ini and [7| allow us to account effectively for both the finite bandwidth of the excitation 
laser and Doppler broadening effects. For this, we will assume excitation by a purely inhomo- 
geneously broadened light source, i.e., the light is described by a monochromatic field with 
fixed irradiance Jq and a statistical probability distribution P{uj) for the frequency. The dis- 
tribution is assumed to be Gaussian, centered at the transition frequency UJ21 and with an rms 
width Aujrms ^ 1/t with r the interaction time, i.e., the time at which the upper state pop- 
ulation is calculated. The Rabi frequency is given by = ^;§^A2i( Ji Mi 1 — g | J2 ^^2)^ -^0 



for a I Ji, Ml) I J2, M2) transition excited by light with polarisation —q = —1, 0, 1. With 
Ji = 1 and J2 = 0, the non-zero squared Clebsch-Gordan coefficients all equal |. Assuming 
the metastable (lower state) atoms to be equally distributed over the three magnetic sub- 
states and the polarisation to be pure but arbitrary, the total upper state population can 
now be evaluated for both the weak and strong excitation limits by integrating Eqs. IHl and 
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P22{t) = —7 — 3-^21^ T weak excitation limit (8) 



(P22) = — i/t — 5—-; strong excitation limit (9) 

6 V n,^21 ^^rms 



BEAM EXPERIMENT 

Here we estimate the feasibility of an in-beam spectroscopic experiment on the 1.557/im 
He* 2^S 1 2 ^Sq forbidden transition. The atomic beam is Zeeman-slowed to 100 ms ^ 
and is transversely cooled to an rms velocity spread of twice the Doppler limit for the He* 
2 ^Si 2 cooling transition {Avrms = 0.6 ms~^). The beam has a diameter of 1 cm and 
an atom flux $ = 10^^ s^^. The transition is excited by the beam of a 2 Watt CW fiber laser 
with a linewidth of 100 kHz, intersecting the atomic beam perpendicularly. For simplicity, 
we will assume a flat-top square intensity profile of size x Dy = 10 x 10 mm. 

In this experiment, the inverse of the interaction time 1/Tint = 10^ s~^, the on-resonance 
Rabi frequency of the transition Q = \j A21 j-^^j-, = 74s~\ and half the decay rate 
r/2 = 25 s^^. We are sufficiently in the weak excitation limit in this case to allow the use of 
Eq. |H| The total upper state population now evaluates to P22(t) = 4.6 x 10^^. This results 
in a flux of excited atoms $e = P22(t)$ = 4.6 x 10^ s~^. 

Using 1083 nm light resonant with the 2^Si 2^P2 transition, after the interaction 
region the non-excited fraction of the atoms can be deflected by simple radiation pressure. 
Using surface ionisation and an electron multiplier on the non-defiected upper state atoms, 
in principle all excited atoms (2 ^Sq) can then be detected. 2 ^Sq atoms produced by the 
discharge beam source can be very efficiently suppressed in the laser cooling stages used to 
prepare the atomic beam. Thus, we can expect a workable signal. 

The spectroscopic linewidth of 400 kHz is dominated by the Doppler width. Decreasing 
the Doppler width as well as reducing the laser linewidth are the ffist steps towards decreas- 
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FIG. 2: Schematic view of proposed setup. 

ing the spectroscopic linewidth. The hmiting homogeneous hnewidth is given by the 10 kHz 
interaction time broadening. Applying multiple crossings of laser- and atomic beam using 
roof-top prisms will not only increase the signal proportionally, but also decrease the inter- 
action time broadening. As a bonus, possible Doppler shifts due to nonorthogonal excitation 
can be monitored and minimised in this configuration. A schematic view of the proposed 
setup, with three crossings depicted, is shown in Fig. |21 

EXPERIMENT IN AN OPTICAL LATTICE 

The 1.557 ;um wavelength of the forbidden transition can also serve well to form a far 
off-resonance optical trap (FORT) for the 2^Si atoms. The dynamic polarisability at this 
wavelength is fully dominated by the contribution of the 2^Si 2^P2 transition that is 
resonant at 1.083 /im. This can be used in a next-generation spectroscopy experiment. As 
the trapping potential is fully insensitive to the wavelength on the scale of a spectroscopic 
scan of the forbidden transition, a single laser can be simultaneously used for trapping 
and spectroscopy. We propose using a cold atomic cloud and transferring this cloud to a 
one-dimensional optical dipole trap formed by a single retroreflected laser beam. 

We assume 5 x 10^ atoms in a cloud with rms size 68 fim (radial) by 340 fim (axial) at 
a temperature of 5 /iK, as typically produced along the route to BEG production |2^. The 
average density in this cloud then equals no = 1.0 x 10^^ cm"^. We now assume a more 
stable 2 Watt fiber laser with a linewidth of 10 kHz. The trapping/spectroscopy laser beam, 
directed along the long axis of the cigar-shaped cloud, has a waist radius of 136 /im. This 
produces a trapping dipole potential with a maximum depth of 24.4 fiK, to which the full 
cloud can be transferred. Analysis of this trap indicates that roughly half the atoms end up 
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in the lowest vibrational state of the standing wave "micro-traps" along the laser beam axis. 
This will effectively cause a strong Doppler-free part in the spectroscopic signal by Lamb- 
Dicke narrowing of the transition. The Rabi frequency = "^^J ^^^21 = 6.28 x 10'^ s~^, 
where the extra factor of two is due to the fact that the atoms are now trapped at the 
antinodes of the standing wave. There is now no a priori fixed interaction time. However, 
as the Rabi frequency is high we can easily choose the excitation to last a time t satisfying 
the time-averaged strong excitation limit {Q ^ 1/t ^ r/2): t ^ 0.2 ms. Eq. Elthen results 
in an averaged excited state fraction (P22) = |\/f a"^^ — ~ '^■^^ ^ 10^^. Given the number of 
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atoms in the trap, every trapped sample will lead to 10^ atoms in the excited state, which 
have to be detected. In an experiment, the trapping/excitation laser can be set off-resonance 
for trapping, and switched to the frequency of the forbidden transition to excite the atoms 
to the 2 ^So state. We can then detect the excited atoms by photoionisation after a few 
milliseconds of excitation (adjusting the geometry of the ionisation laser beam to ionise all 
2 ^So atoms expelled from the trap). 

An alternative detection option is simply by measuring the increased Penning ionisation 
that we expect when atoms are excited. The excited 2 ^Sq atoms can also decay through 
Penning ionisation in collisions with the 2^Si atoms. We assume the rate constant of this 
process to be 10~^° cm'^s"^, i.e., on the order of the rate constant for Penning ionisation of 



unpolarised atoms in the 2 ^Si state 27|. This results in a decay rate of Fp/ = 100 s~^, lead- 
ing to a homogeneous broadening of the transition with this value. However, as the dipole 
shift of the excited 2 ^Sq state has opposite sign, the excited atoms are anti-trapped and 
escape the cloud of 2'^Si atoms in Tesc ~ 0.3 ms. Then Penning ionisation stops, effectively 
reducing the total 2 ^Sq ionisation rate by a factor of 30 as compared to the case of perfect 
overlap. 

As the Doppler width is effectively eliminated by Lamb-Dicke narrowing, the linewidth 
is now determined by the laser linewidth. As this width can be further reduced, ultimately 
the limiting factor will be simply the upper level decay rate. 

However, systematic shifts of the transition frequency will have to be carefully considered 
and corrected for. The largest shift is caused by the difference in dipole shift between 
the lower and upper level of the transition and amounts to 2.0 MHz for the chosen lattice 
parameters. In optical lattice experiments aiming at optical frequency standards one selects 
a lattice laser wavelength at which the Stark shift of the lower state equals the Stark shift 
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of the upper state (the 'magic' wavelength). For hehum there is no practical wavelength 
available. The highest magic wavelengths are around 410 nm, where (accidentally) the dipole 
shift itself is so small that no lattice is feasible at reasonable laser power, and at 351 nm, 
where the polarisability is 15 times smaller than at 1.557 fim and the sign is such that 
the atoms cannot be confined at antinodes. At 1.557 /im, combining measurements of the 
transition frequency at different laser intensities with careful calculations of the intensity- 
dependent dipole shift will still allow for very accurate extrapolation to zero Stark shift. 

Collisional shifts, vanishingly small in the beam experiment, may contribute as well in 
the lattice experiment. For fermions this shift will be absent at the temperatures considered. 



DISCUSSION AND CONCLUSION 



The beam and lattice experiment both promise signal strengths and linewidths that 
should make a measurement with 1 kHz resolution possible. Improvements beyond this level 
also seem feasible. Standard frequency comb technology easily allows an absolute frequency 
measurement at this accuracy. A fiber-laser based frequency comb around 1.5 fim or a 
titanium-sapphire laser based frequency comb (after frequency doubling) may be used for 
this purpose. 

The main obstacles to be overcome are the residual Doppler linewidth for the beam 
experiment and the dipole and collisional shifts for the lattice experiment. Another factor 
to be considered in both experiments is the Zeeman shift due to a stray magnetic field. A 
solution is to measure only anM = O^M = transition. For ^He, this can be achieved by 
simply exciting with linearly polarised light. For '^He, it will be important to shield stray 
magnetic fields and measure both — | ^ — | and +| ^ +| transitions and take the average. 

What will an absolute frequency measurement at or below the 1 kHz level in ^He (or '^He) 

n 

test? In a recent paper Morton, Wu and Drake have tabulated the most up-to-date 
experimental and theoretical ionisation energies, both for ^He and ^He. The theoretical 
uncertainties for the 2^Si and 2^So states are larger than the experimental error for ^He 
for both states. Experiment and theory agree to well within the error bars. It may be 
expected that the theoretical error in the 1.557 fim transition frequency will be smaller than 
the quadratic sum of the errors in the ionisation energy of the metastable states due to 
cancellation effects flsJ]. Therefore, a measurement of the 2^Si - 2 ^Sq transition frequency 
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will not only provide a direct and accurate link between the ortho (triplet) and para (singlet) 
helium system but will also test QED calculations more stringently than the existing data. 

The isotope shift can be calculated with very high accuracy. Using the most recent 
values for the nuclear masses and the most recent evaluation of the difference in the square 
of the nuclear charge radii, i.e., 1.0594(26) fm^ we deduce a theoretical transition 

isotope shift of 8034.3712(9) MHz. A measurement of the isotope shift thus provides a very 
sensitive test of theory. It may also be interpreted as a measurement of the difference in 
nuclear charge radius of '^He and '^He. A measurement at the 1 kHz level will provide this 
difference with an accuracy of 0.001 fm Q|, similar to the accuracy obtained from isotope 
shift measurements on the 2'^Si - 2'^P transition. These have been performed with ~5 kHz 



accuracy 
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16|, limited by the 1.6 MHz natural linewidth of that transition. The 8 Hz 



natural linewidth of the 1.557 fim transition and the possibilities to improve the beam and 
lattice experiment further via sub-Doppler cooling resp. uncoupling of the trapping and 
spectroscopy lasers, may push the accuracy below the 1 kHz level providing new challenges 
to theorists calculating ionisation energies and, in the case of the isotope shift, provide the 
most accurate data on differences in nuclear charge radius and nuclear masses of the isotopes 
involved. 



* Electronic address: Lw,yassen@fe_w,vu,nJ| 

[1] Drake G.W.F. and Martin W.C., Can. J. Phys. 76 (1998) 679. 

[2] Drake G.W.F. , NucL Phys. A737 (2004) 25. 

[3] Eikema K.S.E., Ubachs W., Vassen W. and Hogervorst W., Phys. Rev. A 55 (1997) 1866. 

[4] Bergeson S.D. et al., Phys. Rev. Lett. 80 (1998) 3475. 

[5] Cancio Pastor P. et al., Phys. Rev. Lett. 92 (2004) 023001. 

[6] Pavone F.S. et al., Phys. Rev. Lett. 73 (1994) 42. 

[7] Mueller P. et al., Phys. Rev. Lett. 94 (2005) 133001. 

[8] Sansonetti C.J. and Gillaspy J.D., Phys. Rev. A 45 (1992) Rl. 

[9] Dorrer C. et al., Phys. Rev. Lett. 78 (1997) 3658. 

[10] Lichten W., Shiner D. and Zhou Z-X, Phys. Rev. A 43 (1991) 1663. 

[11] Pachucki K., Phys. Rev. Lett. 97 (2006) 13002. 



10 



[12] Morton D.C., Wu Q. and Drake G.W.F., Can. J. Phys. 84 (2006) 83. 

[13] Pachucki K., Phys. Rev. Lett. 84 (2000) 4561. 

[14] Drake G.W.F., Nortershauser W and Yan Z.-C, Can. J. Phys. 83 (2005) 311. 

[15] Morton D.C., Wu Q. and Drake G.W.F., Phys. Rev. A 73 (2006) 034502. 

[16] Shiner D., Dixson R. and Vedantham V., Phys. Rev. Lett. 74 (1995) 3553. 

[17] Wang L.-B. et al., Phys. Rev. Lett. 93 (2004) 142501. 

[18] Baklanov E.V. and Denisov A.V., Quantum Electronics 27 (1997) 463. 

[19] Pachucki K., private communication, 2006. 

[20] Baldwin K., Contcmp. Phys. 46 (2005) 105. 

[21] Woodworth J.R. and Moos H.W., Phys. Rev. A 12 (1975) 2455. 

[22] Drake G.W.F., Phys. Rev. A 3 (1971) 908. 

[23] Lach G. and Pachucki K., Phys. Rev. A 64 (2001) 042510. 

[24] Lin CD., Johnson W.R and Dalgarno A., Phys. Rev. A 15 (1977) 154. 

[25] Baklanov E.V., Pokasov P.V., Primakov D.Yu. and Denisov A.V., Laser Physics 15 (2005) 
1068. 

[26] Tychkov A.S. et al., Phys. Rev. A 73 (2006) 031603(R). 

[27] Stas R.J.W., McNamara J.M., Hogervorst W. and Vassen W., Phys. Rev. A 73 (2006) 032713. 



11 



